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FIELD OF THE INVENTION 



The present invention relates to a method and to a device for neutralizing, by 
controlled gas injection, liquid slugs or accumulations at the foot of a pipe portion 
greatly inclined to the horizontal or riser connected to a pipe carrying circulating 
multiphase fluids such as hydrocarbons. 



In order to make deep-sea reservoirs or marginal fields sufficiently cost-effective, 
oil companies have to develop new development techniques, as economical as possible. 
It is thus more advantageous to directly transport the two-phase mixture consisting of 
liquid (oil and a little water) and gas in a single pipe or pipeline to onshore facilities in 
order to be separated. A pipe portion greatly inclined to the horizontal (often close to 
the vertical), referred to as riser by specialists, which is connected to the deep-sea pipe, 
is used therefore. However, the gas and the liquid being transported together, flow 
instability phenomena may occur in the zone of connection with the riser, which lead to 
serious development problems. 

In particular, when the gas and liquid inflow rates are low, the liquid phase 
accumulates in the lower parts of the pipeline and stops the gas flowing past. The 
upstream pressure increases and eventually expels the liquid slug to another low part or 
even in the phase separator at the outlet. These accumulation phenomena can reduce the 
productivity and fill pipes designed to receive gas with liquid. One of these phenomena, 
commonly referred to as severe slugging by specialists, has formed the subject of many 
studies, either experimental by means of test loops, or by simulation by means of 
simulation softwares such as, for example, the TACITE simulation code which is 
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notably the object of the following patents or patent applications : US-5,550,761, FR- 
2,756,044 (US-6,028,992) and FR-2,756,045 (US-5,960,187), FR-00/08,200 and FR- 
00/09,889 filed by the applicant. 

This slugging phenomenon is described hereafter in the simple case illustrated in 
5 the accompanying figures where a pipe of low inclination and a riser ended by a 
separator designed to separate the liquid phase from the gas phase are considered. 

The liquid accumulates in the lower part of the pipe and tends to stop the gas 
flowing past. The gas is compressed until the upstream pressure exceeds the pressure 
due to the weight of the accumulated liquid. A long liquid slug is then pushed by the 
4i 10 expanding gas. Under such conditions, an alternating phenomenon is observed, where 
'=f the liquid blocks the gas phase, then flows off under the pressure of the gas and 

H " s eventually accumulates and blocks the gas again. 
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VI More precisely, the periodic process takes place as follows : 
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Stage I : the liquid accumulates at the foot of the riser and stops the gas flowing 
15 past. The pressure rises ; 

Stage II : the upper level of the liquid having reached the top of the riser, the liquid 
phase flows into the separator ; 

Stage III : the gas pocket reaches the foot of the riser and flows into the riser. The 
slug flows into the separator with a much higher velocity ; the gas pocket « explodes » 
20 in the riser ; 
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Stage IV : when the gas pocket reaches the top of the riser, the pressure at the foot 
of the pipe has a minimum value. The liquid falls down along the wall of the riser. It 
accumulates again at the foot of the riser and a new cycle starts. 

A well-known technique referred to as gas lift by specialists allows to overcome 
this phenomenon. It essentially consists in permanently injecting gas at the base of the 
riser to prevent the accumulation of liquid at the bottom of the pipe. Since this 
phenomenon cannot be really controlled, most of the time one is led to inject large 
amounts of gas, which requires considerable compression means. Furthermore, injection 
of large amounts of gas substantially modifies the gas-oil ratio (GOR), which 
complicates the phase separation operations at the top of the riser. 

SUMMARY OF THE INVENTION 

The object of the method according to the invention is to exercise, by modelling the 
instability phenomena described above, an efficient dynamic control over the pressure 
of the gas to be injected into pipes so as to reduce these phenomena as much as 
possible. 

The method according to the invention allows to neutralize, by controlled gas 
injection, the formation of liquid slugs or accumulations at the foot of a pipe portion 
greatly inclined to the horizontal or riser connected to a pipe carrying circulating 
multiphase fluids. This control is essentially exerted by injecting at the base of the riser 
a volume of gas substantially proportional to the mass flow rate variation with time of 
the gas phase of the circulating fluids, and preferably substantially equal thereto, when 
this variation is positive. 
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According to another implementation mode, control is exercised by modulating also 
the volume of gas injected by a quantity proportional to the mass flow rate variation of 
the liquid phase of the circulating fluids, also measured with the same time interval. 

Injection is carried out at any time t, from the mass flow rate variation with time of 
the gas phase of the circulating fluid, measured at a previous time interval 

The implementing device allows to neutralize, by controlled gas injection, the 
formation of liquid slugs or accumulations at the foot of a pipe portion greatly inclined 
to the horizontal or riser connected to a pipe carrying circulating multiphase fluids. It 
comprises gas injection means connected to the base of the riser, means for measuring 
the flow rate of the gas phase of the circulating fluids, and a computer designed to 
control the injection, by the injection means, of a volume of gas substantially 
proportional to the flow rate variation with time of the gas phase of the circulating 
fluids, when this variation is positive. 

The computer is for example suited to control the injection, by the injection means, 
of a volume of gas substantially equal to the flow rate variation with time of the gas 
phase of the circulating fluids. 

According to an embodiment, the device also comprises means for measuring the 
flow rate of the liquid phase circulating in the pipe, the computer being suited to 
modulate the volume of gas injected by a quantity proportional to the measured flow 
rate variation of the liquid phase. 



BRIEF DESCRIPTION OF THE FIGURES 

Figure 1 diagrarnmatically shows a pipeline of low inclination connected to a riser 
greatly inclined to the horizontal, 

- Figure 2 diagrarnmatically shows a steady flow with continuous penetration of gas 
in a gas-liquid separator at the upper end of a riser, 

Figure 3 diagrarnmatically shows the formation of a liquid slug at the foot of the 
riser, 

- Figure 4 shows the stage when the liquid slug reaches the separator at the top of the 
riser, 

Figure 5 diagrarnmatically shows the gas pocket getting into the liquid accumulated 
in the riser, 

- Figure 6 shows the stage when the liquid falls back to the base of the riser, 

- Figure 7 shows the pressure curve in the pipeline at the foot of the riser during the 
previous cycle, 

Figures 8A to 8E diagrarnmatically and respectively show, in a riser, a flow mode 
referred to as bubble flow (8A), a flow mode referred to as intermittent flow (8B), a 
flow mode referred to as churn flow (8C), a flow mode referred to as annular flow (8D) 
and a flow mode referred to as stranded annular flow (8E), 

- Figures 9A to 9G diagrarnmatically and respectively show, in a pipeline, a flow 
mode referred to as stratified flow (9A), a flow mode referred to as wavy stratified flow 
(9B), a flow mode referred to as droplet annular flow (9C), a flow mode referred to as 
dispersed bubble flow (9D), a flow mode referred to as intermittent flow (9E), a flow 



mode referred to as small-pocket flow (9F), and a flow mode referred to as elongated- 
bubble flow (9G), 

Figure 10 shows an example of evolution of the pressure at the foot of a 14-m high 
riser, without any control by gas injection, 

Figures 1 1 A, 1 IB respectively show the flow rate variations of liquid and gas at the 
riser outlet, also without control, 

- Figures 12 and 13 A, 13B correspond to Figures 10 and 1 1 A, 1 IB for a 250-m high 
riser, 

- Figure 14 shows the evolution of the volume fraction of liquid at the riser outlet, 
without control, 

- Figure 15 respectively shows an example of variation of the pressure at the foot of 
the 14-m high riser before gas injection and the stabilization obtained by controlled 
injection according to the method (from 500 s), 

Figures 16A and 16B respectively show, in the same riser, the curves of the 
variation with time of the liquid flow rate and of the gas flow rate, before and after 
control (also from 500 s), 

- Figure 17 shows, in a 250-m high riser, the pressure variation with time at the foot 
of the riser, without control (thick line) and with control (dotted line), 

- Figures 1 8 A, 1 8B respectively show, in the same riser, the liquid and gas pressure 
variations with time, without control (thick line) and with control (dotted line), 

- Figure 19 shows the compared evolutions of the pressure at the foot of a 14-m riser 
with control of the injection pressure as a function of q G and q L (full line) and only of q G 
(dotted line), 
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Figures 20A, 20B respectively show, in the same riser, the liquid and gas flow rate 
variations with time, with control of the injection pressure as a function of q G and q L 
(dotted line) and only of q G (full line), 

- Figure 21 shows, in a 250-m high riser, the pressure variations with time at the foot 
of the riser, with control of the injection pressure as a function of q G and q L (dotted line) 
and only of q G (fall line), 

- Figures 22A, 22B respectively show, in the same riser, the liquid and gas pressure 
variations with time, with control of the injection pressure as a function of q G and q L 
(dotted line) and only of q G (full line), 

- Figures 23A, 23B respectively show the respective liquid and gas flow rate 
variations with time at the riser outlet, and 

Figure 24 diagrammatically shows an embodiment of the device for implementing 
the method, allowing the formation of slugs to be neutralized. 

DETAILED DESCRIPTION 

1-1 Modelling 

Modelling of the flow phenomena in the system consisting of the pipe and of the 
riser of Figure 1 is obtained by means of the following hypotheses. 

We choose a one-dimensional approximation where we consider the averages of the 
various quantities on a (straight) cross-section of the pipeline or of the riser. Since this 
approximation is acceptable only if the radius of curvature of the pipeline is assumed to 
be infinite and its diameter constant, the modelling procedure will concern the parts on 
either side of the connecting elbow. 
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We use a drift flow type modelling procedure with a mass conservation equation 
per phase and wherein the liquid and gas momentum conservation equations are added 
to one another so as to have a single equation referred to as mixture momentum 
conservation equation. To close the system, we choose a relation expressing a friction 
law between the phases. 

We also conventionally assume that : the flow is isothermal, the fluids are 
Newtonian, the gas is a perfect gas, the liquid is incompressible (its density is therefore 
constant) and there is no mass transfer between the two phases. 

We consider that the Mach number of the mixture is small so that the pressure 
waves are propagated at an infinite velocity instead of a velocity close to the sound 
velocity in the mixture. High-frequency phenomena are suppressed but the void fraction 
waves continue to be propagated at a velocity close to the velocity of the mixture. This 
hypothesis is reinforced by the fact that we study the system « in transition » to the state 
of obstruction or disturbed state, i.e. close to the steady state. This hypothesis is 
translated in the model into the absence of inertia terms in the momentum conservation 
equation. 

1-2 Selection of the fundamental quantities 

The following quantities are defined : 

Rg and Rl are the volume fractions of gas and of liquid in the pipes. 
q G and q L are the mass flow rates of gas and of liquid per section unit. 
V G and V L are the velocities of the gas and of the liquid. 
p G and p L are the densities of the gas and of the liquid. 
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P is the mean pressure of the mixture. 

All these quantities are functions of RxR (space and time), zero on RxR_* and 
continuous on RxR,*. According to the context, it is essential to always precisely say if 
the pipe concerned is the pipeline or the riser. Thus, we work on the pipeline in [0,L]xR 
by writing the variables (x,t) and on the riser in [-H,0]xR by writing the variables (z,t). 

However, if an equation is valid in the pipeline as well as in the riser, the equation 
will be formally written with the variables (x,t) e [-H,L] u R+. 

1-3 Intrinsic equations 

All these quantities are connected by algebraic and differential relations which do 
not depend on the flow considered (these flows are described in the next section), which 
is what we call intrinsic equations. The other types of equation, mainly friction laws, are 
studied in the next section. 

1-3-1 Algebraic equations 

We first express the relations directly obtained from the definition of the quantities. 



The perfect gas equation allows to establish the following relation between the 
pressure and the density of the gas : 



<1g =Pg r g v g 



(LI) 



<1l = Pl R l V l 



(1.2) 



R G +R L =1 



(1.3) 



Pg( x >0 = 



a 2 



(1.4) 



In fact, a = I — or Mg is the molar mass of the gas, T the temperature of the 

mixture and R the perfect gas constant, a corresponds to the sound velocity in the gas at 
1 bar. 

1-3-2 Differential equations 

The mass conservation of each phase imposes that : 



^c_ + ^ = 0 (1.5) 
dt dx 

^J^ + ^ = 0 (1-6) 
dt dx 

The momentum conservation according to the previous hypotheses imposes that 



r)P 

?t = -g sm6{R G p G + R L p L )-F p (L7) 

OX 

where 0 is the inclination of the pipe, g the gravity constant and F p the wall friction 
(friction of the stream against the wall). 

A closing equation is added to the aforementioned equations in form of an algebraic 
slippage law as follows : ¥(P, Re, V G , V L , p L ) = 0. 

1-4 Slippage law 

The selected slippage law depends on the flow regime. Three flow types can be 
considered : stratified flow, dispersed bubble (or simply bubble) flow and intermittent 
flow. All the flow regimes are illustrated in Figures 8A to 8E and 9A to 9G. For our 
study, we consider the case of an intermittent flow in the riser and of a stratified flow in 



11 



the pipeline. We shall see hereafter that, considering possible simplifications, no 
slippage law is necessary for the stratified flow regime. 

1-4-1 Intermittent flow 

This flow regime is due to a « superposition » between a bubble flow and a 
stratified flow. When the gas flow rate increases, the bubbles clump together and 
coalesce. Large shell-shaped bubbles appear. They are separated by liquid slugs which 
generally contain small gas bubbles. 

Under such conditions, the friction law for an intermittent flow is expressed as 
follows : 



hence function (V G , V L , V«). 

Besides, V w is experimentally determined and has the following form : 
V m = (O.35sin0 + 0.54 cos 0)y[gD 

II STUDY OF THE PIPELINE-RISER SYSTEM WITHOUT FRICTION 

In this part, we fix F P = 0 in Equation 1.7 considering that the riser is vertical and 
that it is assumed that, in the flow disturbances observed, the frictions only have a 
limited influence in relation to the gravity. Solution of the equations is therefore a priori 
simplified. 



V G -C 0 (R G V o +R L V L )-V a> =0 



(L10) 



II-l Steady state in the riser 

In the case of a steady state, the various quantities only depend on the space ; they 
are marked with a line to show that they are steady quantities. Furthermore, any 
constant function is asterisked. Thus, for example, under steady conditions, the liquid 
and gas flow rates q L and q G are constant. They are therefore denoted by q L * and q G *. 
Besides, the pressure at the top of the riser is denoted by P 0 . 

Intermittent flow 

The closing equation is thus given by L10. Furthermore, as before, in the steady 
state, the gas and liquid flow rates being constant, they are denoted in the same way. 

The implicit formulation of P under intermittent flow conditions is written as 
follows : 



g(P(2)) = g(P 0 )-Z 



(11.15) 



with g given by the relation : 



(l-C 0 )pla 2 q* G 
<1g+ PlK+C 0 q* L 



+ 



SpI(<1g +P*lK + c o?2) 



The pressure can be calculated in implicit form as for a steady flow regime under 
the assumption of an intermittent flow. 



P(z) can thus be assumed to be known at least numerically. It is then posible to 
express all the other quantities as a function of P(z) and z. For the transient states, the 
steady quantities are assumed to be known. 



II-2 Steady model in the pipeline 

We a priori have the flow regime choice. However, practice shows that it is more 
rational to choose the stratified flow regime, or at a pinch the intermittent flow regime, 
knowing that the second one is less probable. The solution of the intermittent flow 

regime is similar to that of the previous section while not imposing {6 = —) in Equation 

1-7. The system can once again be integrated. Assuming that the flow in the pipeline is 
stratified, the closing equation is too complicated to allow the equations to be solved. 
We therefore involve the pipeline only by the effect it has on the riser, i.e. by measuring 
the riser inlet quantities (flow rates, pressure, surface fractions) instead of determining 
them in relation to the calculated pipeline outlet quantities. The calculations are 
therefore replaced by riser inlet measurements. 

II-3 Transient state in the riser in the tangent linear model 

We develop hereafter the calculations about the steady state and linearize the 
intrinsic equations to the first order so as to simply solve the transient system. 

II-3-1 Riemann invariant 

Introduction of this invariant allows to facilitate the pertinent solution of the 
equations in the tangent linear model. It is expressed in the form k = - which 

has no evident physical significance. We thus have the following relation : 



V(z,0 e] - H,0[xR:, ^ + F ™ = 0 (11.28) 

ot oz 



Quantity k moves in the riser, under intermittent flow conditions, at the velocity V G . 

Thus, the property of k allows us to assert that the propagation of k under 
intermittent conditions is expressed as follows : 



Quantity k at a height z of the riser has the same value as at the bottom of the riser 
T(z), one time unit a earlier. 

With this Riemann invariant, we are going to express P'(-H,t), which is the pressure 
variation at the bottom of the riser where the disturbances start, as a function of the riser 
inflow rates. This can allow us to conceive a control over this pressure. 

II-3-2 Transient inflow-outflow laws, intermittent flow regime 

We suppose here that the values of the steady quantities are given (these functions 
are known from the steady mode study). By applying the linearization techniques, we 
can establish the following result : 



k(z 9 t) = k{-H,t-T{z)) 9 



(IL30) 
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(z 9 t) = F(z) N z * q G _ H 



(11.31) 
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where N 3 
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q'c-H , qV-H are functions of t respectively obtained from N, q' G , q' L where the 



first variable has been set at z, -H and -H respectively. 



In fact, linearization of Equation 1.7 gives : 




(IL32) 
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We are going to replace R' G by its expression as a function of P' and k\ We 

therefore note that k = — B^g — g linearizing this relation, we obtain 

(1-C 0 )^ 

= ^(l-Co^)-^/^ if we replace R' G in n.32 by this expression, we 

p G +C 0 k 

obtain the following differential equation : 



dz K a\p G +C 0 k) P G +C 0 k 



The Riemann invariance property of k extends to k' in the following forms 

SL + v G ^- = 0 and k'(z,t) = k'(-H,t-T(z)) with T{z)=[ ^dx 

at oz V q 

We can therefore write : 



10 Furthermore, we can calculate k' as a function of the flow rates, and we obtain 
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By injecting this relation into 11.33, this equation becomes 



dz a 2 (p c +C 0 k) 



g(PG-P*L)(l-CoR G ) k < 



Pc + C 0 k 



<Ig 



kC n 
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q' G (-H,t-T(z)) ^-q\ (-H,t-T(z)) 



We then put : 

I(2) = gRa(C 0 k +P :) ^ = _g(p G -p\)(l-C 0 R G )J_ 
a 2 (p G +C 0 k) P G +C 0 k q G 

This equation is then integrated by means of the constant variation method. We then 
obtain : 
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exp [ H -B(x)dx 

I 



P\z,t) = P\-H,t) + 



kC, 



C(x)\ q' a (-H,t - T(z)) f-q' L {-H,t -T(z)) 

I Pi 



exp[ H B(u)du 



dx 



Finally, we put : 



v) = CoT-\v)el <V ' B ^ du V G oT-\v) and F(z) = e L~ B{x)dx then 



P\z,t) = F(z{ f Z> M(v) q' a (-H,t-v) f-tf L (-H,t-v) 



dv + P\-H,t) 



and if we put N(z,v) = l[ 0 . r(z) ](v)M(v) then 



P'(z,t) = F(z) 



as we have seen above. 



In these relations : 

q 5 G _H is the mass flow rate variation with time of the gas phase in the circulating 
multiphase fluid at the foot of the riser, i.e. at the height -H in relation to the top of the 
riser ; 



q Vh is the mass flow rate variation with time, at the same height, of the liquid phase in 
the multiphase fluid. 

in GAS INJECTION CONTROL 

The action selected to correct the disturbances consists in controlling the pressure at 
the bottom of the riser. In fact, if this pressure remains close to a steady value, this 
means that the slugs do not form and that the gas is never really blocked. The action 
will thus concern the pressure P(-H,t) at the bottom of the riser. 

ni-1 Gas lift control 

We assume that the outlet pressure is fixed and therefore that P 5 0 = 0. Relation 
III. 3 1 can thus be written as follows : 



f Ck ^ 

?G. H — <IL. 
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It can immediately be seen that, if one of the two members of the convolution 

product is zero or very small, it is the same for P'(-H,t). Now, N 0 cannot be modified. 

f 

There still is quantity 
denoted by Q(t). 
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q'G_ H ^-q'L_ H (t) , homogeneous with a flow rate, that is 



The principle of the control mode according to the invention will essentially consist, 
at predetermined action intervals, in measuring at a time t, the above quantity in order to 
obtain a measurement M(tO, then, at the time t, + At, in adding to Q(U + At) the quantity 
u(t, + At) = -M(t,), and so on. We thus have Vi > 1 : 
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+ /AO = -MO, + (/ - l)Af) 
g(f , + fAf) = u(t i + /A/) + + i£t) 

Of course, when At -> 0, it amounts to equating the second term of the convolution 
product to zero. Since it is physically not conceivable to have a control that sucks the 
gas in a two-phase mixture, we use a control such that : 

u(t x + iAt) = max(-M (t x + (i - 1)A/), 0) . 
IV VALIDATION OF GAS INJECTION CONTROLS 

To validate this control mode, we use the aforementioned TACITE code which 
simulates flows in the most accurate and realistic way possible. 

We have studied two very different slug formation instances and tested our control 
mode thereon. 

IV- 1 Comments on simulations 

The TACITE software uses a finite- volume type method to simulate flows in pipes 
and the pipeline is discretized for example according to the gridding method described 
in patent application FR-EN-00/08,200 filed by the applicant. To simulate gas injection 
at the foot of the riser, we modify the flow between the two grid cells situated just 
before and just after the bend, respectively numbered n - 1 and n. The initial state of our 

C k 

simulations is the steady flow regime. Coefficient — — is therefore measured once and 

Pl 

for all at the initial time. We identify the steady flow rates with the pipeline inflow 
rates, and our riser inflow rates with the flow rates of grid cell n - 1. To display the 
results, we approximate to the flow of gas to be injected by the flow rate difference 
between grid cells n and n - 1 . 



IV-2 Simulations of various non-controlled cases 

The simulations carried out by means of the TACITE simulation code are based on 
two geometries where L represents the length of the pipeline, H the height of the riser 
and 0 their common diameter. 

Case 1 :L = 6Om,H=14m,0 = 5 cm, 

Case 2 : L = 1750 m, H = 250 m, 0 = 25 cm. 

Case 1 : 14-m riser 

This is the case which we take as the reference : 

P 0 = 1 bar, Q L = 2.10" 2 kg/s and Q G = 2.10"* kg/s. This case has marked oscillations, 
with an oscillation period of the order of one minute. The flow regime corresponds to 
our hypotheses : stratified in the pipeline and intermittent in the riser. 

The oscillations of the system (Figs. 11 A, 11B) are not great enough to completely 
cancel out the flow of gas at the foot of the riser or for the free surface of the liquid to 
fall below the level of the mouth of the riser. We are therefore not in a situation of 
significant slug formation. This property allows to remain close to our tangent 
linearized hypothesis, while remaining in a flow regime that is not far from the steady 
state. We therefore remain, for any t, in the quasi-steady situation of transition towards 
slug formation. 

Two main stages can be distinguished in the cycle : 

1. A stage of accumulation of the liquid in the riser, the liquid flow rate is zero at the 
outlet and the pressure increases ; 



2. A stage of liquid slug production where the pressure decreases. 
Case 2 : 250-m riser 

The main characteristics of this case are as follows : 
P 0 = 10 bar, Q L = 4 kg/s and Q G = 0.5 kg/s. 

These conditions are closer to the real operating conditions of a pipeline. In the 
absence of any control, the system enters into a liquid slug formation stage. It can be 
observed that the liquid flow rate at the outlet goes through already violent expulsion 
stages, even before the slug expulsion and pressure fall stage. At the end of the 
simulation, the pressure has reached a maximum value, and the system is about to enter 
into the liquid slug expulsion stage. The system is therefore in the phase of transition to 
slug formation only during the first moments, because the starting point of the 
simulation is the steady state. 

Figure 14 shows very fast cycles in the evolution of the liquid fraction. This is a 
sign of high instability in this case. 

IV-3 Control of both cases 

The first control that we are going to study is the theoretical control found in the 
previous chapter. Following our observations, we shall see that it is possible to conceive 
another control, that we are also going to test. 

IV-3-1 Theoretical control test 

We introduce here a control by gas lift by injecting, at the foot of the riser, a flow of 
gas of the following form : 



11(0 = max(0, (<7' G m -^q\ H )(/ - A/)) 

P L 

Case 1 : 14-m riser 

We start control just before t = 500 s. The graphs of FigsiV-16, IV- 17. show in 
parallel the evolution of the free system (thick line) and that of the controlled system 
(fine dotted line). This control allows to maintain the pressure and the outlet flow rates 
close to their steady value. The mass of gas injected is small in relation to what comes 
from the pipeline ; the outlet flow rate increase due to the injection of gas is less than 
5%. 

Case 2 : 250-m riser 

For this simulation, we start control from t = 0 s. The graphs of Figs. 15 and 16 A, 
16B show in parallel the evolution of the free system (thick line) and that of the 
controlled system (fine dotted line). The gas flow rates in the presence of control are 
much more regular than in the previous case. 

Control by gas injection using the theoretical formula thus functions in both cases, 
although they are very different. 

IV-3-2 Simplified control 

Since a great correlation is observed between the evolutions of the gas and liquid 
flow rates, and since we almost always have |#' G a | > 
control the system with the simplified control : 
W (0 = max(0,(4' G „ )('-A0) 



C 0 k 



, it is also possible to 



Case 1 : 14-m riser 

In the graphs of Figs. 17, 18 A, 18B, control with q G and q L is shown with a thick 
line, whereas control with q G only is shown in fine dotted line. Control starts from t = 
Os. 

The results of the control concerning only q G are practically identical to those 
obtained with the control already tested, or even slightly better as regards the pressure 
oscillations and the control speed. 

Case 2 : 250-m riser 

In this distinctly more unstable case, we observe (Figs.22A, 22B) no fast 
oscillations of the gas flow rate with the simplified control. The system controlled with 
q G evolves with a slight lead in relation to the same case controlled with the complete 
expression. 

We can therefore see that the simplified control without coefficient q L also allows to 
control the system, but the theoretical control is better because, unlike the simplified 
control, it does not lead to problems at the riser outlet, it is more economical as regards 
injection gas and the system is controlled just as well. 

IV-3-3 Control robustness 

We test the robustness of our theoretical control (q G and q L ) in relation to the 
reaction time of the sensors and of the actuators. We therefore compare the control 
obtained above by adjusting the gas flow rate at each time interval (thick dotted line) 
with the control obtained by adjusting these parameters with a lower frequency (fine 
dotted line). 
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Case 1 : 14-m riser 

Control is readjusted every 3 seconds here, which corresponds to 20 time intervals 
for TACITE. Control starts from t = 0 s. 

It can be seen in Figs.23 A, 23B that the flow rates are a little more irregular, but the 
system is still controlled. 

Case 2 : 250-m riser 

The system is not stabilized with a 2-second time constant. With a 1 -second time 
constant, the system is stabilized but we remain close to the order of magnitude of the 
calculating interval (about 0.3 s in this case). Furthermore, there is practically no 
difference with the case where control is adjusted at each time interval. 

In any case, we have a good control even though the flow rates oscillate very 
slightly around a mean value. 

The device for implementing the method comprises (Fig.24) gas injection means 1 
connected to the base of the riser, means 2 for measuring the flow rate of the gas phase 
of the circulating fluids, and a computer 3 intended to control injection, through 
injection means 1, of a volume of gas proportional to and preferably substantially close 
to the flow rate variation with time of the gas phase of the circulating fluids, when this 
variation is positive. 



